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Abstract. The Chicxulub impact crater in northwestern Yucatan, Mexico is the primary candidate for 
the proposed impact that caused mass extinctions at the end of the Cretaceous Period. The crater is 
buried by up to a kilometer of Tertiary sediment and the most prominent surface expression is a ring 
of sink holes, known locally as cenotes, mapped with Landsat imagery. This 165 ± 5 km diameter 
Cenote Ring demarcates a boundary between unfractured limestones inside the ring, and fractured 
limestones outside. The boundary forms a barrier to lateral ground water migration, resulting in 
increased flows, dissolution, and collapse thus forming the cenotes. The subsurface geology indicates 
that the fracturing that created the Cenote Ring is related to slumping in the rim of the buried crater, 
differential thicknesses in the rocks overlying the crater, or solution collapse within porous impact 
deposits. The Cenote Ring provides the most accurate position of the Chicxulub crater’s center, and 
the associated faults, fractures, and stratigraphy indicate that the crater may be —240 km in diameter. 


Introduction 

The Chicxulub structure in northern Yucatan, Mexico was first interpreted as an 
impact crater based upon the discovery of a concentric pattern of magnetic and 
gravimetric anomalies (Penfield and Camargo, 1981). A series of recent findings 
indicate that the Chicxulub structure is the long sought after impact crater whose 
iridium-rich ejecta blanketed much of the globe at the Cretaceous/Tertiary (K/T) 
boundary, and which probably caused the mass extinctions that mark this geologic 
transition (Alvarez et al., 1980). Briefly listed, these findings come from studies 
of: (1) mineral grains with shock metamorphism at Chicxulub (Hildebrand el al., 
1991; Sharpton et al . , 1992); (2) circular hydrogeologic anomalies that correlate 
with the magnetic and gravimetric anomalies (Pope el al 1991); (3) the unique 
chemistry of the Chicxulub impact melt rock and its similarity to glassy K/T ejecta 
in Haiti and target rock at Chicxulub (Sigurdsson et al., 1991; Kring and Boynton, 
1992); (4) K/T impact ejecta thicknesses relative to the size and location of 
Chicxulub (Hildebrand and Boynton, 1990; Hildebrand, 1992; Vickery et al., 
1992); (5) impact-generated tsunami wave deposits in the Gulf of Mexico (Hilde- 
brand and Boynton, 1990; Smit et al ., 1992; Alvarez et al . , 1992); and (6) the 
radiometric age (65.0 million years ago) of the Chicxulub impact melt and Haiti 
K/T tektites (Swisher et al 1992). 

In this paper we expand upon our earlier work (Pope et al 1991) and examine 
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Fig. 1. Tectonic features of the Yucatan Peninsula (modified from Weidie, 1985 and interpretation 
of TM imagery). Note locations of exploratory wells (well names given in Figure 6). 


in detail the surface evidence for the buried Chicxulub crater. We combine analyses 
of satellite and geological data to elucidate the unique geology of northern Yuca- 
tan, which we propose developed as a result of an asteroid or comet impact 65 
million years ago. 


Geomorphology 

The northern coastal plain of the Yucatan Peninsula is bounded by the Ticul Fault, 
Gulf of Mexico, and Caribbean Sea (Figure 1). The plain is a region of youthful 
karst development (Isphording, 1975, 1977), and is composed primarily of uplifted 
sections of Early and Late Tertiary platform carbonates, with a band of Quaternary 
beach and lagoon deposits along the coast (Figure 2). There are two main fracture 
systems within the plain: the Holbox Fracture Zone in the east (Weidie, 1982, 
1985), and a system in the west whose salient feature is a semi-circular boundary 
between zones of fractured and unfractured rock (Pope et al 1991) (Figures 1 
and 2). The fracture patterns are interpreted primarily from solution features in 
the Tertiary limestones (Anonymous, 1983). A chain of sink holes is found along 
this boundary zone, named the Cenote Ring (cenote is the local name for sink 
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Fig. 2. Geological map of northwestern Yucatan (Anonymous 1983) showing relationship between 
Cenote Ring and regional fracture patterns. Compiled from sheets F16-7, F15-9-12, and F16-10. 


hole, adapted from the Yucatec Maya word dz'onot (Hall, 1936)). The Cenote 
Ring was first described using Landsat Thematic Mapper (TM) imagery during a 
survey of water resources (Pope and Duller, 1989). 

The Cenote Ring contains over two hundred cenotes, which range in size from 
50 to 500 m wide and 2 to over 120 m deep (Hall, 1936; Soki et al . , 1984). The 
ring is a nearly perfect 165 ± 5 km diameter semi-circle, truncated by the coast, 
and centered on the coast 17 km east of Progreso (Figure 2). Cenote density along 
the ring varies from 3 cenotes/km 2 in a 3 km wide band of cenotes in the southwest 
(Figure 3), to a chain of single cenotes 3 km apart in the southeast. The northeast- 
ern portion of the Cenote Ring is best described as a boundary between a zone 
of abundant cenotes to the east and few to no cenotes to the west (inside) of the 
ring. Southwest of the Cenote Ring is a possible second partial ring. It is composed 
of over 10 cenotes in a concentric arc with a radius of 92 km from the Cenote 
Ring center. 

A topographic profile across the region from the Seasat radar altimeter (Figure 




% 


K. O. POPE ET AL. 



Fig. 3. Landsat TM band 5 (in red), 4 (in green), and 3 (in blue) image (April, 1985) of the 
southwestern portion of the Cenote Ring (location shown in Figure 1). Note chain of black (and 
occasionally blue-green) dots across center of image. These dots are cenotes: collapsed solution cavities 
that intersect the water table. Image width is 15 km. Landsat data from EOSATCo., Lanham, MD, 


4) shows that the three geologic units of the plain form distinct geomorphic 
surfaces. The early Tertiary carbonate surface has a different slope and is much 
more dissected (karstified) than the late Tertiary one. The profile also shows that 
while on the macroscale the late Tertiary surface is the smoothest, it contains a 
series of depressions 5-10 m deep where the profile crosses the Cenote Ring. 
(Figure 4). These depressions contain the cenotes and have a common base level, 
which approximates the regional water table (Marin, 1990). 

Hydrology 

The hydrology of the plain is largely controlled by fractures in the flat Tertiary 
limestones (e.g., Back et al., 1979; Hanshaw and Back, 1980; Back et al., 1986; 
Marin, 1990). The aquifer in northwestern Yucatan is unconfined, except near the 
coast where an impermeable caliche has formed near the surface along the zone 
of fresh water discharge as it migrated inland with rising Holocene sea level (Perry 
et al., 1989). This aquitard varies in width from 5 km near the northeastern 
terminus of the Cenote Ring to 30 km near the southwestern end of the ring. As 
indicated by the fracture pattern noted above, solution features inside the ring are 
rare, but are abundant along and outside of the ring. Except for the partial second 
ring noted above, there are few solution features southwest of the ring. This is 
due possibly to the presence of the impermeable caliche layer, which has sealed 
surface cracks. This paucity of detected fractures may also reflect the fact that 
both rainfall (ca. 500mm/yr) and elevation (less than 8 m above mean sea level) 
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Fig. 4. Topographic profile (location shown in Figure 2) from Seasat radar altimeter. The coastal 
plain is composed of three geomorphic surfaces: a coastal zone with Quaternary beach and lagoon 
sediments, a smooth Late Tertiary surface, and a rough, karstified Early Tertiary surface. The altimeter 
lost lock on the ground surface at the 100 m high escarpment of the Ticul Fault, which demarcates the 
southern boundary of the coastal plain. The altimeter profile falls within the Cenote Ring between 20 
and 70 km inland along track. Seasat altimeter profile is from track No. 43744. Heights are with respect 
to the geoid, therefore elevations above mean sea level are approximate. Seasat data provided by H. 

Frey, NASA Goddard Space Flight Center. 
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Fig. 5. Landsat TM band 6 (thermal infrared) images near Dzilam Bravo, northeastern terminus of 
Cenote Ring (north at top, location shown in Figure 1): (a) wet season, October 1986, (b) dry season, 
April, 1985. Sea-surface cold water anomalies are color-coded black. Light blue areas are warm sea 
water and dark blue areas slightly cooler. Warm land surfaces are in cream, yellow, and red. The gray 
and light blue areas inland from the northern lagoon in the wet season image are clouds and cloud 
shadows. The black anomalies, interpreted as fresh water spring discharge from the Cenote Ring, 
appear within the coastal lagoon and along the seaward side of the beach ridge in the lower left-hand 
(southwest) portion of the wet season (a) image. Spring discharge appears to be greatly reduced in 
the dry season image. The black area offshore is cold, deep, marine water. The TM sensor recorded 
a sea-surface temperature of ~26 °C (light blue), and cold water anomalies (springs) of ~23 °C (black). 
These temperatures are uncorrected for atmospheric effects, and are computed using equations from 
Marlaret et al. (1985). Sea-surface and spring temperatures of 30 °C and 27 °C respectively are reported 
from in situ measurements in the southeastern Yucatan Peninsula (Back et al ., 1979). Differences in 
the absolute temperatures are due to cooling effects of the atmosphere in the satellite measurements, 
however the temperature gradients are similar. Image width is 9 km. Landsat data from EOS AT Co., 

Lanham, MD. 


are the lowest in this region and solution features are greatly subdued. Ground 
water flows toward the coast along these fractures until reaching the ring, which 
acts as a barrier that intercepts the flow. The resulting increased subsurface flow 
along the ring has created a system of caves and conduits whose collapse formed 
the cenotes. 

Ground water flow along the ring is indicated by a valley-shaped depression in 
the ground water surface centered on the ring (Marin, 1990). Additional evidence 
is found along the coast where large fresh water springs occur. Landsat TM thermal 
infrared (band 6) imagery of the northeast terminus of the Cenote Ring near 
Dzilam Bravo contains several cold thermal anomalies along the barrier beach 
and in the brackish lagoon (Figure 5). We interpret these anomalies as discharges 
of cold, fresh water from springs reported from this area (Isphording, 1977). Cold 
water anomalies are very rare along the coast within the Cenote Ring and east of 
the ring terminus. Where they do occur, they are smaller than those near Dzilam 
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Bravo. A major exception occurs at the terminus of the Holbox fracture zone 
where huge cold water anomalies are found, which indicate a high discharge of 
fresh water (Duller and Pope in press). Wet and dry season TM imagery of the 
northeastern ring terminus reveals that spring discharge responds to seasonal rains 
(Figure 5), consistent with water table fluctuations (Marin, 1990). We have not 
analyzed TM imagery of the western ring terminus, but we observed many springs 
there during a 1993 reconnaissance trip. 

Subsurface Geology 

We examined data from nine exploratory oil wells drilled by Petroleos Mexicanos 
(PEMEX) in the vicinity of the Cenote Ring (Figure 1). The well data show that 
the ring interior corresponds with a Late Cretaceous structural low filled with 
undeformed Tertiary sediments (Figure 6). We lack well data from north of the 
ring, although seismic data there do indicate a drop in the K/T surface within the 
suspect crater rim (Hildebrand et al., 1991). Beneath the Tertiary carbonates that 
fill the basin are a series of impact breccias and melt, overlain by probable slump 
deposits. The upper-most Cretaceous strata outside the basin are conglomerates 
and breccias (Weidie, 1985). Underlying this unit is a 2.5 km thick sequence of 
Cretaceous carbonates and evaporites underlain by a thin sequence of Jurrasic 
sediments resting on a Paleozoic metamorphic basement. When the well data are 
plotted as a function of their radial distance from the crater center, offsets in the 
Cretaceous and earlier rocks beneath the conglomerate/breccia unit are apparent. 
These offsets are interpreted as possible normal faults, and one such fault is 
recorded in well Tl. 


Discussion 

Previously, we hypothesized that the Cenote Ring and its associated fracture 
system are the result of post-impact subsidence in the rim of the Chicxulub crater 
(Pope et al ., 1991). The zone of fracturing outside of the ring corresponds with 
our zone of inferred normal faulting, and subsidence along these faults is a possible 
cause of the fractures. These faults may be ring faults with their associated terrace 
zone containing slump blocks. If the Cenote Ring (—170 km dia.) corresponds 
with the crater floor, then a rim diameter of 240 km and a 35-40 km wide terrace 
zone of slump blocks are predicted using the morphometric relationships estab- 
lished for extraterrestrial impacts (Melosh, 1989). These dimensions are consistent 
with our inferred zone of faulting outside the Cenote Ring. 

The inferred fault between the Y1 and Y4 wells, as well as the surface fracturing 
in the same area, are outside our postulated zone of ring faulting. These outer 
features may not be related to the impact structure. Nevertheless, impact craters 
of this magnitude are usually complex structures with multiple concentric features 
(e.g., Croft, 1985; Pike, 1985). Therefore ring faults like those discussed above 
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Fig. 6. Subsurface geology of the Cenote Ring based on well data (well locations in Figure 1). Wells 
are plotted as a function of their radial distance from the center of the Cenote Ring. Stratigraphy and 
description of lithofacies of Yucatan 1 (Yl), Yucatan 2 (Y2), Yucatan 4 (Y4), Yucatan 5A (Y5A), 
and Yucatan 6 (Y6) are from Marshall, 1974 as modified in Weidie, 1985. Lithofacies and correlations 
from Sacapuc 1 (SI), Chicxulub 1 (Cl), and Ticul 1 (Tl) are our own, based on stratigraphic profiles 
presented by Lopez Ramos, 1979, and an unpublished version by G. E. Murray, A. A. Meyerhoff, 
and A. E. Weidie provided to us by A. E. Weidie. Thickness of Tertiary section and information on 
Merida 2 (M2) from Lopez Ramos, 1979). All faults except the one recorded in the Ticul 1 well are 
inferred from offset of beds. Gravimetric profile is from data in Hildebrand et a/., 1991 and is drawn 
through the Cenote Ring center and well Y4. Units 1,2, 3, and 4 were originally designated as Upper 
Cretaceous (Lopez Ramos, 1979; Weidie, 1985), however the radiometric age of 65.0 million years 
(Swisher et al , 1992) for the unit 4 melt rock indicates that units 1 and 2 are K/T boundary or later 
units. Unit 3 contains Upper Cretaceous fossils, but we suggest that this unit is also of impact origin 
and is a K/T boundary unit as well (see text). Key (interpretations in parentheses): 1 - Marl, sandstone, 
shale, and bentonite: the uppermost sediments in Yucatan 6 are globigerinid biomicrites (post-impact 
slump and backwash crater fill deposits). 2 — Bentonitic breccia (impact breccia). 3 - Conglomerate 
and breccia composed of angular to rounded pebble-sized detrital clasts of anhydrite, dolomite, and 
limestone in a dolomitic micrite matrix (reworked continuous deposit from secondary impacts). 4 - 
Andesite, andesitic tuffs, and glass (impact melt and suevite). 5 - Upper Cretaceous mollusc, pelletoid 
micrite and laminated fossiliferous dolomite. The Upper Cretaceous section in the Ticul 1 well contains 
abundant marls and the upper portion may correlate with unit 1 rather than 5.6- Lower Cretaceous 
pelleted rudist biomicrite, biosparite, and milinolid pelmicrite. 7 — Lower Cretaceous laminated dolom- 
ite. 8 - Lower Cretaceous anhydrite. 9 - Lower Cretaceous milinolid dolomite and quartz-silty dolomite. 

10 - Jurassic (?) arenite, shale, anhydrite, and dolomite. 
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are expected and concentric tectonic features at large diameters are also possible. 
The curvature of the Holbox fracture zone (Figure 1) may be partly controlled by 
buried impact structures. 

Given the age of the fractured rocks, the proposed subsidence would have 
occurred in or after the late Tertiary - many millions of years after the impact. 
Long-term tectonic instability of extraterrestrial impact craters has been proposed 
(e.g., Thomas et aL , 1982; Moore and Ahern, 1983), but remains speculative. 
Seismic studies of terrestrial craters indicate that no significant tectonic movement 
occurs after a few million years (Solomon and Duxbury, 1987). These terrestrial 
studies were limited mostly to craters that were much smaller and older than 
Chicxulub, and often deeply eroded, thus the relevance to Chicxulub is uncertain. 
Theoretical models predict that the rim height would reach equilibrium with 
gravity soon after impact (Melosh, 1989). Such models do not preclude subsidence 
in the rim due to burial by Tertiary sediments, or movement along ring faults in 
response to other tectonic forces. 

Although the northern Yucatan Peninsula is not a tectonically active region 
today, the emergence of the northern peninsula and the formation of the Ticul 
fault in the Tertiary indicate that tectonism has occurred during the proposed 
period of crater rim subsidence. Given these regional stresses another mechanism 
that may account for the fracture pattern is differential fracturing due to the 
thickness of the Tertiary rocks. The Tertiary strata overlying the crater floor are 
3 times thicker than the strata outside the rim, therefore these thinner strata are 
much more susceptible to fracturing. The abrupt change in thickness near the 
Cenote Ring may actually concentrate stresses there. 

A non-tectonic hypothesis for the formation of the Cenote Ring was proposed 
by Perry et aL (1992), who suggested that a coral reef may have formed in the 
shallow water of the crater rim, and that solution collapse within the more porous 
reef facies has propagated upward to the surface. If this hypothesis is correct, then 
these early Tertiary reefs must have been much more extensive than Perry et aL 
suggest, because the solution features are not limited to the vicinity of the Cenote 
Ring. The suggestion that the surface solution features are related to the upward 
propagation of deep-seated solution collapse features may, however, be a plausible 
mechanism. Such deep solution collapse could have originated within the highly 
porous conglomerate/breccia unit. We propose that this unit is the product of 
secondary cratering by ejecta, following the process described by Overbeck (1975). 
The final distribution and conglomeritic texture of this unit probably resulted from 
the reworking processes described for marine impacts (Overbeck et al. y 1993). 
The conglomerate/breccia unit is extensive (Figure 6), and its distribution roughly 
corresponds with the surface fracture pattern that created the Cenote Ring. 

Our suggested rim diameter of -240 km is greater than the 180 km proposed 
by Hildebrand et aL (1991), which was based on airborne geophysical data. Infor- 
mation from the Ticul 1 well is critical to the placement of the crater rim; our 
assessment places this well within the crater and Hildebrand and colleagues place 
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it outside. Both the basin morphology and the Late Cretaceous lithologies support 
the placement of the Ticul 1 section within the crater. 

The well data (Figure 6) show that the highest point of the basin rim is close 
to our crater rim estimate of —240 km, and that the K/T boundary in the Ticul 1 
well is depressed about 230 m below this proposed rim. The brief descriptions of 
the Upper Cretaceous strata in Ticul 1 indicate that these strata are unlike Upper 
Cretaceous strata outside the rim (Lopez Ramos, 1979), and apparently the con- 
glomerate/breccia unit in Yucatan 5 A does not extend as far as the Ticul 1 well 
(unfortunately no detailed description of the Ticul 1 well has been published). 
The Upper Cretaceous strata in Ticul 1 include abundant marls, which are not 
present in adjacent wells outside the crater, but are present in the deposits interpre- 
ted as post-impact slump overlying the impact melt (in the Y6, Cl, and SI wells). 
We suggest that the upper portion of the Cretaceous section in Ticul 1 may also 
be post-impact slump deposits. It should be noted that the slump deposits inside 
the crater were originally assigned to the Upper Cretaceous (Lopez Ramos, 1979), 
but radiometric dating of the underlying impact melt (Swisher et al., 1992) com- 
bined with a re-appraisal of the paleontology (Hildebrand et al ., 1991) confirms 
that these sediments are K/T boundary or later deposits. 


Conclusion 

A large body of evidence now exits indicating that a large comet or asteroid impact 
occurred at Chicxulub at the end of the Cretaceous Period. The Cenote Ring and 
its hydrogeological properties add to the unique character of this event in geologi- 
cal history. The hydrogeological data help confirm the circular form of the Chicxu- 
lub structure, and given the nearly perfect circularity of the Cenote Ring and its 
inferred relationship to the buried crater, the ring’s center provides a more accur- 
ate location for the crater’s center than do the geophysical data. 

We identified several possible mechanisms related to the impact that may have 
formed the Cenote Ring. The actual process may well be a combination of these 
mechanisms. Many uncertainties remain about the crater, the most relevant to 
this study are its size and details of its structure. The hydrogeological data shed 
light on these issues and indicate that Chicxulub may be —240 km in diameter, 
which is much larger than previous estimates. Nevertheless, the exact size remains 
uncertain and future research on this well preserved and perhaps largest known 
terrestrial impact crater is needed. 
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